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ABSTRACT 



We examine sequences of dense packings of n congruent non-overlapping disks inside a 
square which follow specific patterns as n increases along certain values, n = n(l), n(2), ...n(k), ... 
Extending and improving previous work of Nurmela and Ostergard [NO] where previous pat- 
terns for n = n{k) of the form /c^, k"^ — 1, A;^ — 3, k{k + 1), and 4A;^ + k were observed, we 
,J^ ' identify new patterns for n = /c^ — 2 and n = k"^ + L^/2j • We also find denser packings than 

C^ . those in [NO] for n =21, 28, 34, 40, 43, 44, 45, and 47. In addition, we produce what we 

conjecture to be optimal packings for n =51, 52, 54, 55, 56, 60, and 61. Finally, for each 
identified sequence n(l),n(2), ...n{k), ... which corresponds to some specific repeated pattern, 
we identify a threshold index fco, for which the packing appears to be optimal for k < ko, but 

^ I for which the packing is not optimal (or does not exist) for k > k^. 

a^ 

en ' 1. Introduction 

\o 

Tjj- ! In a previous paper [GLl], the authors observed the unexpected occurrence of repeating 

"patterns" of dense (and presumably optimal) packings of n equal non-overlapping disks in- 
side an equilateral triangle (see Fig. 1.1 for and example). It is natural to investigate this 
phenomenon for other boundary shapes. In particular, this was done by the authors [LGl] 
for the case of n disks in a circle. However, in contrast to the case of the equilateral triangle 
where the patterns appear to persist for arbitrarily large values of n, for the circle the identified 
packing patterns cease to be optimal as the number of disks exceeds a certain threshold. 

In this note we describe the situation for the square. In a recent paper, Nurmela and 
Ostergard [NO] present various conjectured optimal packings of n equal disks in a square for 
up to 50 disks. They also point out certain patterns that occur there. By using a packing 
procedure different from theirs, we improve on their best packings for n = 21, 28, 34, 40, 43, 
44, 45, and 47. We conjecture that these new packings are optimal, as are the new packings 
we give for n = 51, 52, 54, 55, 56, 60, and 61. 

We confirm all the repeated patterns mentioned in [NO], specifically, for n = /c^, /c^ — 1, 
k'^ — 3, and k{k + l), and we identify two new patterns, namely, for n = k'^ — 2 and n = k'^+[k/2\ 
The latter pattern incorporates the packings of n = Ak'^ + k disks as identified in [NO]. 

It was found in [NO] that the obvious "square" pattern of packings of n = fc^ disks becomes 
non-optimal for n > no = 36. This was done by presenting a configuration of k'^ = 49 disks 




Figure 1.1: The conjectured densest packing of n = 256 disks inside an equilateral triangle, a 
member of the series n = np{k) = A((/c + l)p — 1) + {2p + 1) A(A;) for p = 5 and k = 3, where 
A(m) = m{m + l)/2. The n — p+1 = 252 shaded disks can not move (they are "solid"), the 
p—1 = 4 non-shaded disks are free to move within their local confines (they are "rattlers"). The 
densest packings of n disks for all checked values of the form n = np(k),p = 1, 2, ..., A; = 1, 2, ..., 
have this pattern consisting of one triangle of side {k + l)p — 1 and 2p + 1 triangles of side k 
with p—1 rattlers that are "falling off" the larger triangle. 

with the diameters larger than m = l/{k — 1). The latter m is disk diameter in the packing 
that obeys the pattern. The standard unit of measure used in most papers on the subject is the 
side of the smallest square that contains the centers of disks. We repeat the same procedure 
for the other patterns identified both in [NO] and in the present paper. Namely, for each 
pattern we state the rule of its formation which allows us to compute the corresponding value 
of m = m{n). Then we pinpoint the uq that belongs to the series and such that the packing 
of riQ disks constructed according to the rule is (presumably) optimal, but for which m{ni) 
for the next value rii > uq in the series when the packing is also constructed according to the 
rule is worse than a certain challenger disk configuration (which may or may not be a solid 
packing) . 

In this manner, we confirm observation in [NO] that the best packings oi n = k'^ — 1 disks 
loose their characteristic pattern for n > ni = 48 disks. We also found that although it was 
not stated in [NO], the packing of 47 disks presented there (as well as our better packing of 
47 disks) challenges the series k"^ — 2. Thus, the pattern of the series n = k"^ — 2 becomes 
non-optimal for n > uq = 34, and that of the series n = k'^ — 1 for n > uq = 35. We also 
found challenger disk configurations or packings for other patterns for values of n > 50 which 
were not identified in [NO]. Namely: uq = 56 for the series k{k + 1) (with the challenger 
ni = 72), and uq = 61 for the series /c^ — 3 (with the challenger rii = 78). The situation for the 



series n = n{k) = k"^ + [k/2\ (n = 5, 10, 18, 27, 39, 52, 68,...) is more complex: the pattern 
exists for 5 < n(k) < 52 as a solid packing and is probably optimal for these n(fc) except the 
case n = n(3) = 10 which is the subject of several publications ([G], [Sch], [Schl], [Val]). For 
n = n(8) = 68 the configuration constructed according to the pattern rule has a slight disk 
overlap, i.e., it does not exist as a disk packing, and the overlap persists for all n = n{k) > 68. 
In the geometric packing problem, progress in proving lags that of conjecturing. Thus, 
we should warn the reader that almost all our statements are conjectures; they are based on 
computer experimentations with the so-called "billiards" simulation algorithm [L], [LS]. In all 
series except A;^ — 3, the construction rule we found for generating a pattern for a given n is 
a finite procedure and m{n) can be expressed as the root of a well-defined polynomial. The 
existence of the packing of a given pattern for a fixed n, even if not the pattern's optimality 
(when appropriate), should be considered proven. However, for the series /c^ — 3 we have to 
resort on an infinite simulation procedure [L], [LS] to construct a pattern for an arbitrary n and 
to compute m{n). Hence even the pattern's existence as a solid packing is a conjecture here. 
We point out that some of the proposed methods in the literature attempt to prove a packing 
optimal or, at the least, prove that a packing with particular parameters exists. Usually to 
fulfill this task, the packing must be actually presented, even if only as a conjecture. Thus, 
we try to present our conjectured packings in a verifiable and reproducible form; we provide 
14 decimal digits of accuracy for its parameter m and clearly identify the connectivity pattern 
(touching disk-disk and disk- wall pairs). Some previous papers provided disk coordinates in 
the presented packings. ^ The interested reader can contact the authors directly for the 
coordinates (since that would otherwise take up too much space in the paper). 

2. Packings 

The parameter m supplied with each packing is the ratio of the disk diameter to the side 
length of the smallest square that contains the disk centers. Bonds or contact points mark 
disk-disk or disk-boundary contact. In the packing diagrams, bonds are indicated by black 
dots. Most of the packings presented are conjectures. This means that a proof is required not 
only for their optimality, if any, but even for their existence. Thus, a bond implies a conjecture 
that the corresponding distance is zero while the absence of a bond implies a conjecture that 
the distance is strictly positive. We placed or did not place a bond between two disks or a disk 
and a boundary based on the numerical evidence: the bond was placed when the corresponding 
distance was less than 10~^^ of the disk diameter. Such a choice of a threshold is supported 
by the existence of a well-formed gap between a bond and a no-bond situation: In all cases 
when the bond is not present between apparently touching surfaces, the computed distance 
is at least 10^^ of the disk diameter, and, except for the packing of 47 disks in Fig. 2. 6, it is 
at least 10^^. The existence of this gap also testifies to that in all the packings the double 



^We were unable to reproduce the best packing of 21 disks for which [MFP] provides the diameter m = 
0.27181675 but no other data, e.g. no contact diagram. Our best packing of 21 disks has a smaller diameter 
(see Fig. 2.1). 



precision resolution we employed for the computations sufficed. All solid disks, i.e., those that 
can not move, are shaded in the packing diagrams; the non-shaded disks are rattlers — they are 
free to move within their confines. Different shadings of disks in some packings and a unique 
numerical label for each disk on a diagram are provided to facilitate the discussion. These are 
not part of the packings. 




21 disks 
m = 0.27181225535931 39 bonds 
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23 disks 
m = 0.25881904510252 56 bonds 
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22 disks 
m = 0.26795840155072 43 bonds 
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24 disks (2,2) 
m = 0.25433309503025 56 bonds 




24 disks (2,3) 
m = 0.25433309503025 56 bonds 



24 disks (3,3) 
m = 0.25433309503025 56 bonds 



Figure 2.1: The densest packings found of 21 to 24 disks 




25 disks 

m = 0.25000000000000 60 bonds 




27 disks 

m = 0.23584952830142 55 bonds 
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26 disks 

m = 0.23873475724122 56 bonds 
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28 disks 

m = 0.23053549364267 57 bonds 




29 disks 

m = 0.22688290074421 65 bonds 



30 disks 

m = 0.22450296453 109 65 bonds 



Figure 2.2: The densest packings found of 25 to 30 disks 




31 disks 

m = 0.21754729161912 55 bonds 
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33 disks 

m = 0.21 132838414326 65 bonds 
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34 disks (2,5;2,4) 

m = 0.20560464675957 80 bonds 
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32 disks 

m = 0.21308235294443 61 bonds 
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34 disks (2,4;2,4) 

m = 0.20560464675957 80 bonds 
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34 disks (3,5;2,4) 

m = 0.20560464675957 80 bonds 



Figure 2.3: The densest packings found of 31 to 34 disks 
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34 disks (2,5;2,5) 

m= 0.20560464675957 80 bonds 
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35 disks (2,3) 

m= 0.20276360086323 80 bonds 
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35 disks (2,2) 

m= 0.20276360086323 80 bonds 
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35 disks (3,3) 

m= 0.20276360086323 80 bonds 
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36 disks 

m= 0.20000000000000 84 bonds 



37 disks 

m= 0.19623810145141 73 bonds 



Figure 2.4: The densest packings found of 34 to 37 disks 
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38 disks 

m= 0.19534230412691 77 bonds 





,^ v,^ v,^ v.^ 

39 disks 

m= 0.19436506316151 80 bonds 
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40 disks 

m= 0.18817552201832 85 bonds 



41 disks 

m= 0.18609951 184812 100 bonds 




42 disks 

m= 0.1842770721 1710 90 bonds 
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43 disks 

m= 0. 1 80 1 9 1 1 3545743 85 bonds 



Figure 2.5: The densest packings found of 38 to 43 disks 
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44 disks 

m = 0.17863924567120 82 bonds 
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45 disks 

m = 0.17571631417559 94 bonds 




46 disks 

m = 0.17445936087241 91 bonds 




47 disks 

m = 0.17126830721141 94 bonds 




48 disks 

m = 0.16938210954876 101 bonds 




49 disks 

m = 0.16738607686833 120 bonds 



Figure 2.6: The densest packings found of 44 to 49 disks 
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50 disks 

m = 0.16645462588286 104 bonds 
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51 disks 

m = 0.16561837431260 99 bonds 




52 disks 
m = 0.16538623796964 105 bonds 




54 disks 

m = 0.15913951630719 115 bonds 




55 disks 

m = 0. 1 5755574752972 1 1 3 bonds 




56 disks 

m = 0.15615650046215 119 bonds 



Figure 2.7: The densest packings found of 50 to 52, and 54 to 56 disks 
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60 disks 

m = 0.14950565404867 129 bonds 




61 disks 

m = 0. 1 48544 1 26695 18 121 bonds 




62 disks (3,7 ;2,5) 

m = 0.14569394327531 140 bonds 
not the best packing 




72 disks 

m = 0.13549029317569 169 bonds 

not the best packing 




72 disks 

m = 0.13541666666667 152 bonds 

not the best packing 




Figure 2.8: The densest packings found of 60 and 61 disks and inferior packings of 62, 72, and 
78 disks 
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5 disks 
m = 0.70710678118655 12 bonds 




10 disks 

m = 0.4 1 666666666667 2 1 bonds 
not the best packing 




10 disks 

m = 0.42127954398390 21 bonds 




13 disks 

m = 0.36609600769643 25 bonds 




18 disks 

m = 0.30046260628867 38 bonds 
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78 disks 

m = 0.13046077259640 
not a rigid packing 



Figure 2.9: The densest packings of 5, 10, 13, and 18 disks, an inferior packing of 10 disks, 
and a challenger configuration of 78 disks 



14 
3. Series of packings with similar patterns 

n = k^ . The packings of this series have an obvious square pattern with disk diameter 
m = l/{k — 1). The pattern is optimal for n = 1,4,9, and 16, and probably 25 and 36. For 
ki = 7 this yields value of tti = 1/6 which is smaller than the experimental diameter of the 
challenger packing of 49 disks in Fig. 2.6. 

n = k^ — 1 . The pattern can be viewed as a square of (A; — 1)^ disks arranged in straight 
rows and columns into which one row and one column of "shifted" disks are inserted. A 
"shifted" row or column each has k — 1 disks. In the diagrams for fc = 5 (n = 24, Fig. 2.1) 
and /c = 6 (n = 35, Fig. 2.4), the "straight" disks are shaded more heavily than the "shifted" 
ones. Several equivalent packings are obtained by inserting the shifted row and column in 
different places among the straight ones. A pair of insertion position (i,j),2 < i,j < k — 1, 
identifies the packing, e.g., in packing of 24 disks (2,3) (Fig. 2.1), the shifted row is the second 
and the shifted column is the third (counting from the top left corner). Among the (A; — 2)^ 
packings thus obtained, there are many congruent pairs. The pattern is fully developed at 
n = 8 (one packing) and remains optimal for n = 15 (one packing), n = 24 (three equivalent 
packings), and n = no = 35 (three equivalent packings). For rii = 48 the pattern looses its 
optimality, which can be shown as follows. The angle at the center of disk 13 in packing (2,2) 
of 24 disks (see Fig. 2.1 ) in the triangle formed by centers of disks 13, 19 and 3 is 15° with 



cos 15° = .5\/2 + v3. Similar angles are the same in all the packings of the series which easily 



implies m = l/{k — 3 + y 2 + \/3)- For ki = 7 this yields value of m = 0.168581424 which is 
smaller than the experimental diameter of the challenger packing of 48 disks in Fig. 2.6. 

n = k^ — 2 . This pattern is similar to the pattern of series k'^ — 1, only here there are two 
shifted rows and two shifted columns. The (non-optimal) packing of 62 disks depicted in Fig. 
2.8, shows this pattern for k = 8. There are (k — 3){k — 4)/2 possible ways to insert a pair of 
shifted rows (11,12) because of the restrictions 1 < ii < 12 < k and similarly for the columns 
(ji, J2). Hence, there are ((/c — 3)(/c — 4)/2)^ possible different index sets (ii, i2',ji,J2) for each 
of which we can construct an equivalent packing. Many of these are congruent. The pattern 
is fully developed at k = 5 (n = 23, one packing) and remains optimal for only one more 



value k = kf) = 6 (no = 34, four equivalent packings.). Here m = l/{k — 5 + 2\/ 2 + v3). For 
k = ki = 7 this yields value m = 0.1705406887 which is smaller than that of the challenger 
packing of rii = 47 disks in Fig. 2.6. 

n = k^ — 3 . This pattern is represented by conjectured dense packings of 22, 33, 46, and 
61 disks (Figs. 2.1, 2.3, 2.6, 2.8). A feature of the pattern is a (shaded more heavily on 
the pictures) densely packed "straight" square of (k — 3)^ disks in the bottom left corner. It 
follows that the pattern is not optimal for A; > 10 (n > 97) because the "straight" square 
itself isn't. However, even for a smaller n = ni = 78 we were able to get a disk configuration 
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that challenges the pattern as presented in Fig. 2.8, although this configuration is a not fully 
formed solid disk packing (see Fig. 2.9). 

Looking at the (non-optimal) packing of 78 disks (see Fig. 2.8), the pattern can be further 
described as three alternating columns at the right and three alternating rows at the top with 
one rattler at the top right corner. In each of the three additional rows and the three additional 
columns, most disks are not touching each other. The exceptions are pairs 42-49 in the bottom 
additional row, 73-77 in the top additional row and 76-72 in the right additional column (Fig. 
2.8). Also there is almost full contact for pairs of disks between adjacent additional rows 
and columns, except for the pair 42-75 which are not touching between the first and second 
additional rows. These features are identical in all the packings of the pattern. 

Each such packing can be obtained by "tightening" a certain configuration C described as 
follows: In (7, disks 77 and 76 are not touching each other but all disks in additional rows 
and columns are touching their neighbors. {C is not a solid packing.) We take C as an initial 
condition for the "billiards" packing algorithm [L], [LS]. To "tighten" the configuration, the 
"billiards" algorithm allows the disks to move chaotically in the square without overlaps while 
their diameter increases at a common rate until no further growth is possible. It is remarkable, 
that for each n = /c^ — 3, k = 5, 6, 7, 8 and 9, this chaotic negotiation always converges 
(independently of the initial velocities) to the described pattern with parameter m identical to 
double precision in all runs involving the same n. Moreover, when instead of the configuration 
C, we start with a random initial configuration and zero initial diameters of the disks, the 
same configuration with the same precision results in the runs that achieve the largest m, - 
but only for /c = 5, 6, 7 and 8, excluding value /c = 9. Also, for /c = 4 (n = 13) which seems 
to be the smallest n for which the pattern may exist, both methods, the one beginning with a 
zero-diameter random configuration and the one beginning with the configuration C, lead to 
the same packing which, however, deviates somewhat from the described pattern, see Fig. 2.9 

n = k(k + 1) . The pattern consists of A; + 1 alternating columns with k disks each, see, e.g., 
the packing of 30 disks in Fig. 2.2. Following the path of disks 2, 12, 6, 18, 3, 29 in this 
packing, we come to the expression k cos a for the horizontal side of the square that contains 
the disk centers, where the angle a = a{k) is that at the vertex 2 in the triangle 12, 2, 6. The 
vertical side of the square is A; — 1 + sin a, where the latter expression results from the path 12, 
2, 5, 30, 24, 7. Equating these two expressions for the side we determine from the resulting 
equation: cosq;(A;) = (fc^ — k + \/2k)/{k'^ + 1) and m = {k cos a{k))^^ . 

Looking again at the packing of 30 disks in Fig. 2.2, circles 1 and 2 are separated by a 
non-negative distance because a{k) < 30° or cos a{k) > \/3/2 for k = 5. The latter inequality 
holds and the pattern exists as a solid packing of non overlapping disks only for A; > 4, i.e., for 
n = 20, 30, 42, 56, 72,.... 

We show that the packing of the pattern is not optimal for n > ni = 72 as follows. We 
present another regular pattern of packings of k{k + 1) disks; for k = 8 this alternative pattern 
is depicted in the second row and second column in Fig. 2.8 (with one rattler). The alternative 



16 

pattern exists for all values of n for which the main pattern exists. The diameter m = fh{k) 
of a disk in the alternative pattern is given by m{k) = l/{{k + 1/2) cos f3{k) + ( v3/2) sin [3{k)) 
where (in the example in Fig. 2.8) (3{k) for /c = 8 is the angle at disk 56 in the triangle 67, 56, 
42. As before, we equate the horizontal and vertical sides of the square: k cos /3+cos(/3+7r/3) = 
{k — 1) sin(/3 + vr/3) + sin/3. From this equation we easily determine (3 = P{k) and then we 
find that fh > m for all n > rii = 72 (but fh < m for n = 20, 30, 42, and 56). Note that the 
alternative packing, although it is better than the packing of the pattern for n = 72 disks, is 
not optimal. For example, for 72 disks we found experimentally a disk configuration with an 
irregular structure (not shown) that is better than both the main or the alternative pattern. 

n = k^ + [k/2j . The pattern, as exemplified by cases k = 5 {n = 27, Fig. 2.2) and k = 6 
(n = 39, Fig. 2.5), consists of A; + 1 alternating columns, odd columns having k disks each and 
even columns having k — 1 disks each. As before, with a denoting the angle, say, at disk 15 in 
the triangle 7, 15, 10 in the packing of 27 disks in Fig. 2.2, we compute the side-length of the 
square obtained in two different ways: A; cos a = 2{k — l)sina. From this equation we easily 
determine cos a{k) and then m = 1 / (k cos a{k)) . The pattern exists when, as in the example 
of 27 disks (Fig. 2.2), disks 12 and 13 are separated by a non-negative distance. This occurs 
when sin a(/c) > 1/2. Thus, the pattern exists only for A; = 2, 3, 4, 5, 6, and 7 (n = 5, 10, 18, 
27, 39, and 52) and those packings (shown in Figs 2.9, 2.2, 2.5, and 2.8) are, indeed, optimal, 
except the case of n = 10: proved [GMPW] for n = 5, 10, 18 and conjectured for n = 27, 39, 
and 52. 

For k > 8 the "ideal" pattern yields overlaps in pairs of disks. The overlap increases with 
k. For A; = 8 (n = 68) the overlap is less than 1% of the disk diameter so it is naturally to 
expect that the optimal packing will be a small deviation from the main pattern. Experiments, 
indeed show that best configurations are of this sort. Unfortunately, none of them is a solid 
packing, because it becomes very difficult to find exactly in which pairs the disks are touching, 
and not merely just very close to each other. 
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